Topology induced anomalous defect production by crossing a quantum critical point 
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We study the influence of topology on the quench dynamics of a system driven across a quantum 
critical point. We show how the appearance of certain edge states, which fully characterise the 
topology of the system, dramatically modifies the process of defect production during the crossing 
of the critical point. Interestingly enough, the density of defects is no longer described by the Kibble- 
Zurek scaling, but determined instead by the non-universal topological features of the system. Edge 
states are shown to be robust against defect production, which highlights their topological nature. 
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Phase transitions occur between two different states of 
matter with different orders. When a critical point is 
crossed with a certain rate, the order can be established 
only partially, and the system displays a non-vanishing 
density of defects. The Kibble-Zurek (KZ) mechanism is 
the paradigm to analyze the dynamics of a system driven 
across a phase transition [l|, Q. Accordingly, the density 
of defects produced during the crossing is uniquely de- 
termined by the universality class of the system. Re- 
markably enough, this mechanism accurately describes 
also quantum phase transitions, occurring at absolute 
zero temperature @, S, H, i, 0, B E K EH; the effect of 
crossing a quantum critical line was shown in T^, 13| . 
Here we show that the topology of the system may in- 
duce an anomalous defect production which strongly de- 
viates from the KZ scaling law. We demonstrate that 
such deviation is caused by the properties of the charac- 
teristic edge states, which are a fingerprint of the system 
non-trivial topology. 

The effect of boundary conditions is peculiar for the 
class of systems possessing a topological order [T3]- In 
fact, such systems are characterized by a ground state 
degeneracy that strongly depends on the topology of the 
of the manifold on which they live. Besides, a hallmark 
of topological order is the appearance of states localised 
at the boundary of the system, the so-called edge states. 
Edge states emerge naturally in a wide variety of sys - 
tcms, such as the fractional quantum Hall effect [1 51 ] - 
one-dimensional spin models [16| . or Majorana fermions 
on a lattice [lj} • Recently, edge states have been realized 
experimentally in topolo gica l insulators connected with 
quantum spin Hall effect HI, 19| ■ 

The subject of this letter is to investigate the effects 
of topology in the out-of-cquilibrium dynamics induced 
by the crossing of a quantum critical point. In this con- 
text, the preparation of topological order via adiabatic 
evolution was addressed in [2(|. Moreover, the scaling 
of defects created by crossing a line of quantum critical 
points of topological nature was studied in [l3| ■ The sce- 
nario that emerged supports the KZ mechanism, with a 



generalized scaling law that takes into account the infi- 
nite number of critical points crossed. Let us remark that 
no peculiar behavior induced by the topological nature of 
the system has been found so far. In the following we will 
show how edge states dramatically modify the KZ scaling 
law, thus indicating the need of a new paradigm to de- 
scribe the dynamics of topological defect production. To 
address this problem we consider a specific model orig- 



inally proposed by Creutz in [2l|. As we shall discuss 
below, this system offers a rich scenario where to study 
the influence of topology on the defect production across 
a quantum critical point, and may serve as a paradig- 
matic model for more general topological systems. 

The Creutz model describes the dynamics of a spinlcss 
electron moving in the ladder system depicted in Fig. [TJ 
as dictated by the following Hamiltonian 



H := 



r ie a\ l+1 a n +c l0 b\ l+1 b n 



(1) 



+K (b] l+1 a n + al +1 b n ) + Ma f n b n + h.c 



where we have introduced the fcrmionic operators a n , 
and b n associated to the n—th site of the upper and 
lower chain respectively. The hopping along horizon- 
tal and diagonal links is described by the parameter K, 
and the vertical one by M; additionally a magnetic flux 
9 G [— 7r/2, 7t/2] is induced by a perpendicular magnetic 
field. 

For small vertical hopping M < 2K, the model has a 
second order quantum phase transition at 6 C = with 
equilibrium critical exponents v = z = 1 (sec Eq. [3]). 
The quantum critical point separates two band insula- 
tors with topologically distinct configuration of currents 
(leading to opposite circulation of thereof along the lad- 
der). We thus focus on a quench of the magnetic flux 



6(t) 



jt — 7r/2 G [— 7r/2,7r/2] at constant rate v q . Ac- 



cording to the KZ mechanism, as the system crosses the 
critical point and the energy gap vanishes, excitations are 



unavoidably produced as 



KZ 



CX Vg 



where d is the 



2 



system dimension. We now explore the influence of the 
ladder topology on the dynamics across a quantum crit- 
ical point and describe the circumstances where strong 
deviations from the KZ scaling arise. 

Let us first discuss the quench dynamics of a ladder 
with periodic boundary conditions cll+i = cti, = b\, 
which impose a closed topology on the system. In- 
troducing the fcrmionic operators in momentum space 



y 



One" 



and b„ 



where q 



— 7 

T, J 



with j = 1, L, the Hamiltonian in equation JT]) reads 
i? = -2if^( a t,fct) 



cos(q — 6) m + cos q 
m + cos q cos(q + 9) 



(2) 

where m := M/2K is the relative vertical hopping, and 
K fixes the energy scale. Hence, the phase diagram de- 
pends upon the parameters m and 9 with the following 
energy spectrum E q ± := E/2K 



E q ± = — cos q cos 9 ± 



sm 2 gsin 2 i 



+ (m + cosq) 2 . (3) 



In the following we will focus on the regime with m = 
(sec fig. 2(a)). Let us notice however that the results 



obtained are robust under small fluctuations of m < 1. 

We shall consider an initial ground state correspond- 
ing to 9(0) = —ir/2 and localised within a plaquette at 



(n,n + l), | - 2K), 



,t 



bl 



l n+l 



4+i)|0). 



Such state arises due to the destructive interference be- 
tween different hopping paths, which forbids electron 
tunneling to next-nearest neighboring rungs n n ± 2. 
At the end of the quench 6(tf) = +7r/2, the system 
shall be found in a state ^(ff)) given by the superpo- 
sition of the ground state still localised within the pla- 
quette (adiabatic evolution of the initial state), and de- 
localized defects which have been excited close to the 
critical point. This dynamics, schematically depicted in 
fig. 3(a) can be rigorously described as a collection of 
L Landau-Zcner processes [22|, [23| (see eq. ©). The 
total density of excitations at the end of the quench is 

Pdei = E^oK-^l*^))! 2 ' (^1 bcin S thc eigenvectors of 
the Hamiltonian. 

In fig. 4(a) we present the density of produced defects 
as the thermodynamic limit is approached. It readily 




FIG. 1: Diagram of the Creutz ladder with periodic bound- 
ary conditions: Electrons freely move across the ladder hop- 
ing along horizontal, vertical and diagonal links that join the 
different sites of the lattice. 



follows that the density of defects is accurately described 
by the scaling P^ef oc which is in complete agree- 

ment with the KZ prediction for the critical exponents of 
the model. In fact, despite the non linear nature of the 
quench of 9, equation© can be linearized close to gap- 
less mode q g = tt/2 at 9 C , yielding thus the conventional 
scenario where the KZ scaling holds. 

Once the periodic ladder has been thoroughly de- 
scribed, we may study the influence of a different topol- 
ogy on the quench dynamics of the ladder. In the case of 
open boundary conditions a^ + i = 0, = 0, thc trans- 
lational invariance of the system is lost. The dynamics 
can no longer be described as a collection of L uncoupled 
two-level systems associated to a given mode q, and tran- 
sitions between all energy levels occur. We observe that 



the quantum phase transitions in fig. 2(b) still resembles 
the periodic case in fig. 2(a) However, the opening of the 
chain is crucially accompanied by the appearance of two 
new levels whose energies lie close to zero. Such levels, 
which correspond to edge states pinned at the ends of the 
lattice, are a direct consequence of thc ladder topology 
(i.e. whether it is an open quasi-one-dimcnsional chain, 
or a closed quasi-one-dimcnsional ring). Below, we study 
the effects of these new states, and thus the effects of 
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(a) Energy spectrum for the periodic ladder 
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(b) Energy spectrum for the open ladder 

FIG. 2: Energy spectrum of the Creutz ladder for m = as 
a function of the magnetic flux 6: (a) For periodic boundary 
conditions we observe that the energy gap vanishes at 6 C = 
0, which correspond to a quantum critical point, (b) For 
open boundary conditions we have two additional zero-energy 
levels which correspond to edge states that characterize the 
topology of the ladder. 
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(a) Defect production in the periodic ladder 




(b) Defect production in the open ladder 

FIG. 3: Scheme of the production of delocalised defects in the 
magnetic flux quench 6(t): (a) In the periodic case, starting 
from the localised plaquette | — 2K) n , delocalised defects are 
produced as the system is driven across the critical point 9 C = 
0. At the end of the quench 6 — > tt/2, the ladder will be in 
a linear superposition of the plaquette-like ground state and 
de-localised excited states, (b) For open boundary conditions, 
we have to consider the dynamical quench of the topological 
edge state localised at the left boundary of the ladder. 



topology, in the ladder quench dynamics. 

For open boundary conditions, we can initialise the 
system in two qualitatively different localised states, 
namely, the plaquette-like state localised at the bulk 
(or boundary), or the topological edge state located at 
either side of the ladder \l) := ^(a>\ — i&i)|0), and 



|r) := 7^(-i4 + &L)|0) (sec fig. [3(b)] ). In the case in 
which the system is initialised in a plaquette-like state, 
the density of defects scales according to the KZ mecha- 
nism Pdof c< -JVq- Being plaquette-like states indepen- 
dent on the ladder topology, we consistently find the 
same scaling obtained in the periodic case. Moreover, KZ 
scaling holds regardless of the plaquette position within 
the ladder, either located at the bulk (fig. |4(b)[ ) or edge 



(fig. 4(c)). These results rule out the existence of a 



non-topological attraction mechanism at the boundaries, 
which would modify the ladder quench dynamics. 

Let us now focus on the quench dynamics of a topolog- 
ical edge state localised in the left side of the ladder |Z). 
Interestingly enough, the energy gap between this edge 
state and the remaining excitations shows a similar scal- 
ing as the gap of a localised plaquette, and thus the KZ 
mechanism would results in a similar scaling Pjgf oc ^/v^- 
Contrarily, we have found (see Fig. |4(d)[ ) the surprising 
result that the scaling of the density of defects for an 
initial topological edge state deviates pronouncedly from 
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(a) PBC: Dynamical quench of a bulk plaquette 
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(b) OBC: Dynamical quench of a bulk plaquette 
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(c) OBC: Dynamical quench of an edge plaquette 
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(d) OBC: Dynamical quench of an edge state 

FIG. 4: Scaling of the density of defects produced during a 
magnetic flux quench 9(t). Periodic boundary conditions: (a) 
density of defects Pdef produced after the quench of the mag- 
netic flux as a function of the quench rate v q in logarithmic 
scale for a ladder with L G {80, 160, 320, 640, 1280} rungs for 
an initial plaquette-like state. In the thermodynamic limit, 
the density of defects fulfills the KZ scaling Pd c f oc y/Vq~. 
Open boundary conditions: (b) plaquette localised in the bulk 
\-2K) L/2 for L e {10,20,40,80, 160} rungs, (c) plaquette lo- 
calised in the boundaries | - 2K)i for L G {10, 20, 40, 80, 160} 
rungs, (d) edge state \l) for L G {10, 20, 40, 80, 160} rungs . 



4 



the KZ theory 



-Pdef oc v, 



1.35 



(4) 
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Accordingly, edge states are more robust (with respect 
to non-topological states) because a smaller amount of 
defects is generated during adiabatic quenches c g < 1. 
This anomalous scaling is due to topological constraints 
that decouple the edge states from the excitations that 
lie close to the gapless mode, whose universal density of 
states gives rise to the KZ scaling. Therefore, we can 
state that the quench dynamics of these topological edge 
states is not dictated by the universality class of the sys- 
tem, but rather induced by non-universal effects of the 
topology. It is precisely this fact which does not allow 
a theoretical prediction of Eq. ((4| based on the critical 
properties of the quantum phase transition. 

Summarizing, in this letter we have presented a de- 
tailed study of the effect of the topology on the quench 
dynamics. We have shown that the influence of edge 
states, which constitute a clear signature of the system 
topology, leads to a peculiar scaling of quench defects. In 
this respect, we have presented an insightful model, the 
so-called Creutz ladder, where the influence of these edge 
states on the quench dynamics can be neatly described. 
As we noticed, an adiabatic quench of the magnetic flux 
produces a smaller amount of defects when the ladder is 
initially in a topological edge state. Being the presence 
of in-gap edge states a generic feature of topological band 



insulators [2J, |25j, |26| , we expect that an anomaly in the 
defect production could emerge beyond the specific model 
we considered. We would like to conclude the letter com- 
menting on the relevance of our results for quantum in- 
formation. The suppression of the defects production 
law suggests that the edge states are particularly suitable 
candidates to implement a protected quantum memory, 
even in the presence of strong fluctuation of a magnetic 
field that may drive the system across a critical point. 
Let us finally note that a natural candidate to experi- 
mentally implement this topological quench dynamics is 
that of spinor Fermi gases in optical lattices [23]. The 
two species of fermionic operators in the ladder may cor- 
respond to the fermionic operators associated to different 
spin-components, whereas the non-interacting regime is 
reached by means of Feschbach resonances. 
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